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1.  Introduction 


Various  theoretical  scattering  models  have  been  developed  over  the  years  for 
characterizing  the  electromagnetic  responses  of  trees  within  the  context  of 
microwave  remote  sensing  and  foliage-penetrating  radar  applications.1^15  These 
research  efforts  can  be  loosely  categorized  into  2  main  classes:  The  first  class  is 
derived  from  radiative  transfer  theory,  which  assumes  the  tree  canopy  to  be  a 
uniform  layer  containing  a  random  distribution  of  scatterers  and  relies  on  the  use 
of  statistical  averages  for  the  calculation  of  extinction,  source,  and  phase  matrices 
(or  coefficients);  the  second  class  is  based  on  the  so-called  discrete  scatterer 
approach,  which  is  a  more  deterministic  method  that  seeks  to  directly  account  for 
individual  scattering  mechanisms.  (Note  that  a  vegetation  canopy  can  also  be 
treated  as  a  continuous  medium  with  a  fluctuating  permittivity  function;  this 
approach,  however,  is  most  appropriate  and  tractable  for  weakly  scattering  media 
in  which  the  fluctuating  function  is  small  relative  to  its  mean  value.5)  While 
traditional  radiative  transfer  theory  has  been  applied  successfully  to  a  number  of 
remote  sensing  scenarios  involving  complex  forest  canopy  structures,16  this 
method  does  not  account  for  the  coherent  scattering  effects  occurring  within  the 
tree  structure  and  therefore  cannot  be  used  to  compute  the  phase  response  of  a 
scene.  On  the  contrary,  the  discrete  scatterer  approach — which  first  decomposes  a 
tree  structure  (trunk,  branch  complex,  and  leaves/needles)  into  simple-shaped 
scatterers  such  as  finite-length  dielectric  cylinders  and  dielectric  disks  and  then 
calculates  the  collective  return  by  summing  the  scattered  fields  from  each 
scatterer — is  a  coherent  technique  that  can  be  employed  to  deduce  the  generalized 
polarimetric  response.  The  development  of  the  discrete  scatterer  approach,  then,  is 
dependent  upon  a  thorough  understanding  of  the  scattering  effects  of  dielectric 
cylinders  and  disks — preferably  through  the  derivation  of  closed- form  expressions 
for  their  scattering  matrices.  The  subject  of  interest  in  this  work  is  the  scattering 
from  a  finite-length,  dielectric  circular  cylinder,  which  is  the  elementary 
component  used  to  model  the  trunk  and  branching  structure  of  a  tree.  (The 
leaves/needles  of  a  tree  are  commonly  approximated  as  disks.  At  the  lower 
frequency  bands,  the  effects  of  the  leaves/needles  are  often  assumed  to  be 
negligible;  thus,  a  study  of  the  scattering  from  disks  is  not  included  in  the  current 
work.) 

The  organization  of  the  work  is  as  follows.  In  Section  2  a  complete  derivation  of 
the  scattering  matrix  of  the  finite-length,  dielectric  circular  cylinder  is  presented. 
A  closed-form,  approximate  solution  is  obtained  by  applying  the  modal  (or  eigen¬ 
function)  expansion  method  in  conjunction  with  the  volumetric  equivalence 
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principle;  the  formulation  exploits  the  supposition  that  the  internal  fields  of  the 
finite-length  cylinder  are  the  same  as  those  of  the  infinite-length  case.  The  theory 
is  first  developed  for  a  cylinder  in  free  space.  Subsequently,  the  problem  of  the 
cylinder  located  above  a  finite-conducting  ground  is  considered  in  Section  3,  in 
which  a  multiray  technique  is  employed  to  deal  with  the  effects  of  the  half  space; 
essentially,  this  approach  calculates  the  total  scattering  response  by  coherently 
summing  the  first-order  direct  response  and  the  single  and  double  ground  bounce 
contributions.  An  arbitrarily  oriented  cylinder  is  treated  by  augmenting  the 
solution  with  appropriate  coordinate  transformation  matrices.  In  Section  4, 
simulation  results  are  presented  to  demonstrate  the  utility  of  the  analytical 
solution.  Finally,  in  Section  5,  a  summary  of  the  overall  work  is  given. 


2.  Semi-Exact  Solution  in  Free  Space 


In  the  current  problem,  the  incident  and  scattered  wave  directions  ( kt ,  ks )  and 
their  associated  polarization  vectors  ( h  ,  v)  are  defined  by 


kt  =  sin  6i  cos  +  sin  0j  sin  q \y  -  cos  0tz  ; 
ks  =  sin  Gs  cos  (j):x  +  sin  0s  sin  (j)sy  +  cos  0sz ; 


Ks  = 


zxk . 


z  x  k , 


v.  -h.  xk.  ; 

l,S  l,S  l,S  ? 


(1) 

(2) 

(3) 

(4) 


in  which  the  subscripts  i  and  s  identify  quantities  related  to  the  incident  and 
scattered  waves,  respectively;  0.  is  the  incidence  angle  in  elevation  (measured 
from  -  z  );  6S  is  the  scattering  angle  in  elevation  (measured  from  +  z  );  and  $  y  is 

the  incidence/scattering  angle  in  azimuth  (measured  from  +Jc).  Without  the  loss 
of  generality,  a  vertically  polarized  ( TMZ )  plane  wave  is  taken  as  the  source  of 
excitation;  that  is,  with  time  convention  of  ej(0t  and  designating  the  wave  number 
in  free  space  as  ka ,  the  incident  electric  field  is  simply  E'  (r )  =  E ovie~jk°kt'1  .  This 

wave  is  impinging  on  a  dielectric  circular  cylinder  with  radius  a  and  relative 
dielectric  constant  er.  Here,  for  the  initial  discussions,  the  cylinder  is  centered  at 
the  origin,  with  its  axis  aligned  with  the  z-axis.  The  first  step  is  to  derive  the  fields 
inside  an  infinite-length  cylinder — which  has  an  exact  modal  solution; 
subsequently,  the  scattering  response  of  the  finite-length  case  is  deduced  by 
approximating  the  cylinder’s  internal  fields  with  those  of  the  infinite-length  case. 
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It  can  be  shown  that  the  z-component  of  the  electric  field  of  the  incident  wave  can 
be  expanded  as 

+oo 

K  (  A  4>,z)  =  Ea  X  f"Jn  ( kpip )  eMi) ,  (5) 

n=— oo 

with  Ea  =  -Eo  sin  6ejk:,z ,  k_  .  =  ka  cos  0, ,  and  kp ;  =  ka  sin  9i .  Note  that  //'  =  0  for 

a  TMZ  wave.  It  is  theorized  that  the  z-components  of  the  scattered  and  internal 
electric  and  magnetic  fields  of  the  infinite-length  cylinder  can  be  written  in  the 
following  forms: 


k(p,*,z)=e.  £  ; 

n=—co 

(6) 

+00 

E‘(p,lz)  =  Ec  Y,  B„J,(K.rP)e  • 

n=— oo 

(7) 

+CO 

//:  (p,  !#,  z)  =  E.  Y  C,H?  (kMp) e1^ ; 

«=— oo 

(8) 

+00 

Hi  (p.lz)  =  E„Y  D,J,  (*,,rp)+MI 

(9) 

n=— oo 


where  k  =  k0yjsr  -  cos2  ;  and  ,  Bn,  Cn,  and  Dn  are  unknowns  to  be 
determined. 


Noting  that 


1  — 


H, 


1  — 


1 

8Ez 

1 

oH  _ 

(10) 

or  pa  ^ 

jco2  pa 

P 

d(f> 

j(os 

5/7 

i  2  N 

K,i 

-k  . 

Z,l 

1 

dHs 

1 

.  _L 

dEz 

(11) 

or  pa 

jm1  pa 

P 

d(/) 

| 

h 

dp 

the  ^-components  of  the  same  set  of  electric  and  magnetic  fields  can  be 
formulated  as  follows: 


K(p,^,z)  =  -e0 


cos  6: 


Kp  sin-  Oi  t 


-t-ou 

X  j~nnJn(kp,ip)e 


(12) 
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E;(p,iz)=Eofj 


Edt(p,</>,z)  =  E„2 


-  cos  6j 
k.p-za  0, 


”-i,H(;'(kp,p)- 


jas„  sin  4 


p) 


Jn(<l-<ti). 


(13) 


-cos  <9, 


koP(£r~coS2  &;) 


nBnJn(kp,rP)- 


j(OS0yjsr -cos2  e, 


DJn'{kP,rP) 


o. 


H‘t{p,t,z)  =  E. 

+00 

h;(p,  *,+=£„£ 


j(D/u0  sin  0i  „= 
+00  1  —  cos  0, 


+00 

Y,rj,iK,P)e 


jn(0-<k)  . 


(14) 

(15) 


kaP  sin- 3 


jcaii0  sin  0j 


a,h?'  iKip) 


(16) 


Hi,(p^.z)  =  E0fj 


-COS  <9, 


koP(^r  -  COS2  0,.) 


nD„Jn(kP,rp)- 


sk„ 


JOJli,  yjcr  -cos' (l 


bj:(K.,p) 


Jn(  (H4) 


(17) 


The  unknowns  An ,  Bn ,  Cn ,  and  Dn  can  be  found  by  matching  the  tangential 
electric  and  magnetic-field  components  at  the  cylinder’s  surface.  It  can  be  shown 
that 


where 


A.,  =■ 


H 


(2) 


{kP.ia) 


—M 

R 


~Jn(kp,ia) 


-r+lM„  . 
’  jMp^k 


B,  = 


C  =D 

n  n 


j,{kP,a)  . 


-j  "k0nW cos 6* 

OMoJn(kP,ra)Rn  ’ 


(18) 


(19) 


(20) 


(21) 
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R 


7i (k  a)  H(2)  (k  .a)  r 

=  1  P'  ’  1  P'  1  •  [M,  •  Nn  -  W2n2  cos2  0.  ] ; 


Mn=Hn{kpia)-Jn[kpra)- 

Nn=Hn(kpia)-£jn(kpa)- 


{kpja)  = 


H?\kpJa) 

kp,aH^(kp,aY 


j  tk  a\-  jApA  . 


w  = 


( kpJa )~  ( kpra ) 


(22) 

(23) 

(24) 

(25) 

(26) 

(27) 


Consequently,  the  p  -component  of  the  electric  field  inside  the  cylinder  can  be 
expressed  as 


Edp(P’<P’z)  =  EoY, 


'  nDn  " 

K(OSrS0P; 


f 

1 

V 


cos"  6, 


s..  -  cos"  0, 


Jn{kp,rp)- 


ij 


jB„  cos  0, 
y]e,-cos 2  ei 


Jn  '(kP,rP 


(28) 


The  fields  inside  an  infinite-length  cylinder  Ed(p,<j>,z )  are  then  completely 
characterized  by  Eqs.  7,  14,  and  28. 


In  view  of  the  volumetric  equivalence  principle,  an  approximate  solution  for  the 
scattering  response  Es  (r,0s,<f>s)  of  a  cylinder  with  length  L  can  be  calculated 

using  the  internal  fields  of  the  infinite-length  problem  treated  above;  that  is,  in  the 
far  field, 


j  | 


5 


—jk0r  +LH  a  2 n 


Anr 


-1/2  0  0 


(29) 


After  substituting  Eqs.  7,  14,  and  28  into  29,  the  integration  over  the  cylinder’s 
volume  can  be  facilitated  by  noting  the  following  relations: 


+1/2 


JVV'(  cos0i+cos^z,  =  Zsinc 


kaL{  cos  6 i  +  cos  6S )  \ 


(30) 


-1/2 


cos 


'  =  -2  nf 


j  cos  </>sJn  '{kp  sp ')  +  nsm&  ') 


Kp,P 


I  -  r  ■ 


(31) 


J”{^.  )  _  FT  . 

C  °/7’ 


(32) 


1 


4;r 


n/d 


KppP 


s- 


na 


KkP-rP 


Jn{kp,rP')-jb'Jn'{kprp') 


c. 


(33) 


i  a 

1  r 


4;r 


J 


7lZ>' 


i“'j.iK,P'y— J,(kpy) 

Kp.rP 


C„  + 


na 


KkP.rP 


-,J„(K.rP')-AJ,:(KrP') 


=  /[(*'->')  ]  A*  -*  > ; 


(34) 


(35) 


where 


p'  dp 


p'  dp 


6 


7”  =  T - ¥~Ti - V  (kpsaJ»  (^aK+i  (kp,ra)~ kp,saJ„  (kp,ra)Jn+ 1  (**,«)) 


;  (36) 


kp,s=koSin0s- 


(37) 


After  much  algebraic  manipulation,  the  scattering  matrix  elements  corresponding 
to  a  vertically  polarized  excitation  wave  are  shown  to  be  given  by 


S'pv  (U)  =  Ps  - -  smc 


ko(£r-l)L  ■  ( k0k (cos 0:  +  cos 0s ) 


V 


+00 

S[k,,+A.+A.. 


(38) 


where  p  =  v,h  and 


(39) 

K2,  =  -lkJ-f  -  (P,  +  ja„)(x+jy)  ■ 

Zlip,r 

(40) 

K3,„  =  i„rj ; 

(41) 

_  jMn  sin  Oi  cos  6t  . 

"  jAkAr. 

(42) 

a  nW  sin6(cos6( 

■f.(K,a)R.  ’ 

(43) 

_  .iMn  sin  ei 
'  )R.' 

(44) 

For  a  horizontally  polarized  excitation  wave,  the  derivation  for  the  scattering 
matrix  elements  S'ph  [ks,k^  closely  follows  the  one  outlined  above;  it  is  seen  that 
they  can  be  obtained  from  Eq.  38  simply  by  redefining  an ,  /3n ,  and  yn  as 
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a, 


(45) 


-nWsin8i cos  0j 

Jn(kp,ra)Rn 


Pn 


jN„  sin  Oj  , 

Jn(kp,ra)Rn  ’ 


(46) 


-n  W  sin  6i  cos  8i 

Jn(kp,ra)Rn 


(47) 


Expressions  for  and  S‘ph  \  ks,k  \  here  match  the  formulations  from 

Karam  et  al.17  (Note  that  there  are  typographical  errors  in  Eq.  25  of  that  journal 
article.17) 


3.  Multiray  Solution  in  the  Presence  of  a  Half  Space 


For  the  case  when  the  cylinder  is  located  above  a  finite-conducting  ground  (with 
the  air-ground  interface  at  z  =  0),  a  total  scattering  matrix  can  be  constructed  by 
summing  the  direct  wave  from  the  cylinder  and  the  ground-reflected  waves. 
Succinctly,  the  first  4  main  contributions  can  be  written  as 


S(U)-F  +7"  (U,)-F(*,KA(i^ 


+r 


(*.)'F  M) 


e-jK{kr?o-Krrg)  +r 


(48) 


where  fo  is  the  location  of  the  center  of  the  cylinder;  r  =  ro  -2{ro  -z)z ; 
kt  g  =  ki  - 2 ikj  ■z\z\ks g  =  ks  - 2 (ks  ■  z ) z ;  S  (•,•)  is  the  cylinder  response  in  free 
space;  and  T^kjs  j  is  the  ground-reflection-coefficient  matrix 


r 


(49) 


with  rv\ki  I  and  rh  ykt  s  j  as  the  ordinary  horizontally  and  vertically  polarized 

ground  reflection  coefficients,  respectively.  The  first  term  in  Eq.  48  is  the  direct 
wave  from  the  cylinder;  the  second  and  third  terms  correspond  to  the  single 
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ground  bounce  returns;  and  the  last  term  is  the  double  ground  bounce 
contribution. 


Note  that  if  the  global  coordinate  system  is  aligned  with  the  local  coordinate 
system  used  for  the  definition  of  S  (•,•)  in  Section  2,  then  S  (•,•)  =  S  (•,•). 

=S  . 

Otherwise,  for  an  arbitrarily  oriented  cylinder,  the  scattering  matrix  S  (•,•) 
should  be  modified  as 


=g 

5 


/\  AO l 

V  -Vs 

v  -hgl' 

1 

1 

II 

_ 1 

1 

1 

’  ^ 

i  i 

*  * 
hf-h,  \ 

(50) 


in  which  ( v.  ,  hj  ),  ( v, ,  hs  ),  and  (  ka ,  kb )  are  the  original  global  coordinate 
vectors;  ( vf ,  hf)  and  (  vf ,  hf  )  are  the  polarization  vectors  in  the  cylinder’s 
local  coordinate  system  written  in  global  coordinates — these  correspond  to  k'a 

*  1  *  ,  ’=^8  /V  *  J  =4>  /V 

and  kh ,  respectively.  Note  that  ka  =  T  ka  and  kb  =  T  kh .  Here,  it  is  assumed  that 
vector  transformations  from  global  to  local  coordinates  (or  vice  versa)  can  be 

=ig  =gi  f =ig  \T 

accomplished  with  a  transformation  matrix  T  (or  T  =  T  ). 


4.  Numerical  Results 


The  backscattering  cross  section  of  a  tilted,  finite-length,  dielectric  circular 
cylinder  is  shown  in  Figs.  1-3  for  the  free-space  and  half-space  cases,  as 
calculated  with  Eqs.  38  and  48.  The  cylinder  has  radius  a  =  15  cm  and  length 
L  =  1.5  m  and  is  centered  at  fo  « (0.49,1.26,3.58) ,  with  its  axis  parallel  to  the 

vector  (0.12,0.32,0.94).  The  cylinder  and  ground  have  relative  dielectric 
constant  and  conductivity  (er,  ad)  of  (13.9,  39  mS/m)  and  (5.45,  20  mS/m), 


respectively.  The  backscattering  cross  section  is  defined  as  a  -  An 


S.J-, 


(~K,k) 


The  incidence  angles  here  are  set  as  0 \  =  40°  and  (j>.  =  200° .  Because  of  the 


dihedral-like  effect,  an  enhancement  in  the  return  is  observed  for  the  co-polarized 
responses  when  the  cylinder  is  located  above  a  ground;  a  similar  effect  is  not 
apparent  for  the  cross-polarized  responses  in  this  example.  Note  that  reciprocity 
dictates  Shv  =  -Svh  in  the  backscattering  direction.18  Although  the  formulation 


given  in  Section  2  obeys  this  relation,  the  multiray  approach  given  in  Section  3 
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does  not.  In  view  of  this  deficiency  in  the  half-space  solution,  the  average  cross- 
polarized  response  (that  is,  (Shv  -  Svh  )  is  plotted  in  Fig.  3 . 


Fig.  1  Backscattering  cross  section  (as  function  of  frequency)  of  a  tilted,  finite-length, 
dielectric  circular  cylinder  located  in  free  space  and  above  a  half  space:  vv  response. 
Parameters:  a  =  15  cm;  L  =  7.5  m;  0.  =  40° ;  and  <f>.  =  200°  • 


Frequency  [MHz] 

Fig.  2  Backscattering  cross  section  (as  function  of  frequency)  of  a  tilted,  finite-length, 
dielectric  circular  cylinder  located  in  free  space  and  above  a  half  space:  hit  response. 
Parameters  are  the  same  as  those  in  Fig.  1. 
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Fig.  3  Backscattering  cross  section  (as  function  of  frequency)  of  a  tilted,  finite-length, 
dielectric  circular  cylinder  located  in  free  space  and  above  a  half  space:  hv  response. 
Parameters  are  the  same  as  those  in  Fig.  1. 

5.  Conclusions 


The  complete  derivation  of  a  closed-form  formulation  for  the  scattering  matrix  of 
the  finite-length,  dielectric  circular  cylinder  has  been  presented.  The  semi-exact 
solution  is  obtained  by  approximating  the  internal  fields  of  the  finite-length 
cylinder  with  those  of  the  infinite-length  case,  which  has  an  exact  solution 
available  through  modal  analysis.  Application  of  the  volumetric  equivalence 
principle  enables  the  generalized  far-field  response  of  the  scatterer  to  be 
characterized  in  an  efficient  manner.  The  theory  is  first  developed  for  a  cylinder 
in  free  space;  subsequently,  the  problem  of  the  cylinder  located  above  a  flat 
dielectric  ground  is  considered  with  a  multiray  technique  in  which  the  ground 
effects  are  taken  into  account  through  the  introduction  of  reflection  coefficient 
matrices  and  phase  delay  terms.  The  scattering  response  of  an  arbitrarily  oriented 
cylinder  is  obtained  by  supplementing  the  canonical  formulation  with  appropriate 
coordinate  transformation  matrices.  The  complete  solution  established  above 
facilitates  the  development  of  the  discrete  scatterer  approach  for  characterizing 
the  scattering  return  from  tree  structures.  An  evaluation  of  the  accuracy  of  the 
solution  and  a  determination  of  its  region  of  validity  is  the  subject  of  Part  II  of  this 
study. 
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